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Abstract: We present a new model for the dynamics of the presynaptic intracellular 
calcium concentration in neurons evoked by various stimulation protocols. The aim 
of the model is twofold: We want to discuss the calcium transients during and after 
specific stimulation protocols as they are used to induce long-term-depression and 
long-term-potentiation. In addition we would like to provide a general tool which 
allows the comparison of different calcium experiments. This may help to draw 
conclusions on a wider base in future. 

Introduction 

A most fascinating challenge in neuron biology is a deep understanding of the mecha- 
nisms involved in long-term-effects (LTE) such as long-term-potentiation (LTP) and 
long-term-depression (LTD). The multiplicity of possibly important mechanisms is 
immense and is explored in a great variety of experiments |IJ. However, the interpre- 
tations of those experiments are not as conclusive as they could be if it was possible 
to better compare experiments executed on different systems. We claim that a lot 
of detailed information on LTE is hidden in presently available experiments. 

One possible way to uncover this hidden knowledge is to construct a tool which is 
able to translate different experiments into a common language and then to compare 
them quantitatively. Such a tool is provided here for the analysis of the intracellular 
calcium concentration in presynaptic nerve terminals. It is very well established by 
experiments that LTE are mostly connected with dynamical changes in the calcium 
concentration H. Larger calcium concentrations induce a lot of mechanisms that 
may influence the transmission efficiency of synapses. This applies not only to the 
postsynaptic side of the synapse, where an calcium influx is believed to be necessary 
for the induction of LTP ||. Intracellular calcium also induces exocytosis on the 
presynaptic side of the synapse and therefore is in parts responsible for the release of 
neurotransmitter into the synaptic cleft M. It is believed that calcium/calmodulin 
dependent protein kinases may trigger the amount of released neurotransmitter. 
This has been shown for example for the calcium/calmodulin kinase II in the squid 
giant synapse ||. Therefore, we think that the intracellular calcium concentration 
is an appropriate observable to study the induction of effects as LTD and LTP from 
the very beginning. 

Already on this early level of LTE-induction a comparative quantitative evalua- 
tion of different experiments may lead to new insights. Corresponding experiments 
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have been done using different systems or using the same system under different 
conditions. In addition, a dependence of the presynaptic calcium dynamics on the 
corresponding target cell has been observed ||. In order to compare those exper- 
iments quantitatively one has to determine characteristics of the experiment that 
are important for the calcium dynamics and to introduce them into the model ter- 
minology. 

In the following we will develop a corresponding model which is based on pre- 
viously developed models of the intracellular calcium dynamics. The number of 
models describing presynaptic calcium dynamics is comparably small |7], & £| |TU|. 
The main focus of model work in this field lies on the postsynaptic calcium dy- 
namics [0, |TTJ, [12], [13], [14], [15|, |16|, [17], |l8l (for a review we refer to [19|]). This is 



surely related to the fact that dendrites are better accessible in experiment than 
presynaptic boutons. Therefore, we use the knowledge which has been established 
on the postsynaptic side of the synapse in order to build a new model which can 
be used for comparison of presynaptic calcium dynamics in different experiments. 
Presynaptic measurements became more frequent due to technical developments in 



the last years [20 



The model has to be adjusted to specific experiments in a well-defined proce- 
dure, which is illustrated for the example of presynaptic nerve terminals in the rat 
neocortex 0. The results found here are discussed in the context of LTE. For more 
details of the model and the results presented here we refer to [^TJ. This especially 
concerns a more detailed analysis of the model assumptions and the robustness of 
the results. 



The model on the level of single proteins 

We construct a new deterministic one-compartment model for the presynaptic cal- 
cium dynamics. The intracellular calcium concentration is changed in response to 
variations of the membrane potential by an ion-flux through voltage-gated calcium 
channels. The flux is driven by the electro-chemical gradient at the membrane. In- 
tracellular calcium is bound by an endogenous buffer. Only a small number of ions 
remains free and increases the concentration of free calcium. An increased calcium 
concentration is reduced by two types of transport proteins: ATP-driven calcium 
pumps and natrium-calcium exchangers. Their activity basically depends on the 
free calcium concentration. A leakage transmembrane calcium current is assumed, 
that compensates the activity of transport proteins and channels in the rest state 
of the cell. 

The model is formulated in terms of a set of ordinary coupled differential equa- 
tions for the intracellular calcium concentration c: 

dc G PvJu(gv(U),U(t),U(c)) - p P Jp(g P (c)) - p E J E (g E (c))+L 

It ~ If i + e 6 (c) + 6i(c) (1) 

where Pu,p,e are the surface densities of the voltage-gated channels (HVA), the 
PMCA-type calcium pumps, and the natrium-calcium exchanger (type 1), respec- 
tively. 

>Aj,P,E = gU.P.E -fu,P,E (2) 
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are the corresponding single protein currents. They are determined by the the single 
protein open probabilities <7u,p,e, where gp^ depend on c according to standard Hill 
equations ]22]| for the pumps and the exchangers with Hill coefficient 2 [23] and 



1 p4| , respectively. The open probability for the voltage gated calcium channels 
(HVA) obeys [T2§ 

n -l 



dgy(U) 1 
dt t 



1 + exp 



U 1/2 - U 



9v(U) , (3) 



with t = lms the channel mean open time, U1/2 = 3mV the half activation voltage, 
and k = 8mV the steepness factor. The values are taken from a representative 



HVA-channel ||25|| . I\j is the voltage dependent current which is assumed to follow 
Ohms law with open channel conductivity of YlpS. The electrochemical gradient is 
calculated with respect to U(c), the (calcium dependent) calcium reversal potential, 
which depends on c according to the Nernst equation. 

Ip = 10~ 17 A and Ie = 40 10~ 17 A are the maximum activities of the transport 
proteins. G is a geometry factor (ratio of surface and volume in the compartment), 
z = 2 is the valence of the calcium ions, F the Faraday constant, and L the leak 
current which is determined by the steady state condition. 

U(t) is the stimulating transmembrane voltage function. We assume that the 
feedback effect of the calcium influx on the membrane potential is small with re- 
spect to standard action potentials. One may think of a voltage clamped situation. 
However, for large calcium concentrations as they may appear in high frequency 
stimulations, calcium buffers may saturate and the impact on the membrane po- 
tential may become non- negligible. Such effects are not considered in the present 
analysis. 

The buffers (endogenous and indicator) are treated in a quasi-steady state ap- 
proximation, which claims that the calcium binds and dissociates faster than the 
typical time scale under consideration. Then the dynamical behavior of the buffers 
reduces to a correction factor in Eq. ([!]) which depends on the calcium concentration 
only: 

[K h + c) (Ki + c) 

Here 6 max and i max are the total concentrations of the endogenous buffer and the 
indicator, respectively. K bji are the corresponding dissociation constants. Note that 
it is important to include the indicator used in the experiment, because the indicator 
is basically acting as additional buffer and may drastically change the amount of 
free intracellular calcium. 

The differential equations are solved numerically and the solution describes the 
time course of the calcium concentration resulting in response to single action po- 
tentials or to series of action potentials (as they are used to induce LTD or LTP). 
The above mentioned aspired generality of the model is reflected in a separation of 
the model parameters into three classes, described in the following subsections. 

Universality 

The model is universal enough to be applicable to a wide class of different neuron 
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types. To this end the model is based on the experimental knowledge about single 
proteins which are postulated to have neuron-type independent properties. These 
are the single protein characteristics g x (including all derived physiological prop- 
erties) and the single protein transmembrane currents J x , where x stands for the 
corresponding type of protein (see e.g. G3, 123, 



Type specificity 

The model is specific enough to be applicable to well defined neuron types. This is 
achieved by the introduction of measurable neuron-type specific parameters which 
has to be determined for each experiment separately. Basically, these are the protein 
densities p x in the membrane. As no space resolution of the calcium concentration is 
considered, these densities may be thought as average values over the whole synaptic 
membrane. Also the concentration of the endogenous buffer 6 max and its dissociation 
constant Kf, belong to the neuron-type specific class of parameters. Finally, the 
surface to volume ratio G of the synaptic compartment quantitatively determine 
the concentration changes due to transmembrane currents. 



Condition specificity 

The model includes enough general specifications in order to adjust the model to 
specific experimental conditions. The form and amplitude of the action potential 
U(t) is simulated with a system of coupled differential equations (not shown here) 
and can be adapted to the specific action potential used in experiment. The LTE- 
stimulation protocols used in experiment are simulated with a corresponding series of 
single action potentials. Intracellular calcium concentrations are generally visualized 
with the help of calcium indicators. They act as an additional buffer in the cell and, 
therefore, may influence the calcium dynamics. In the model they are treated in 
complete analogy to the endogenous buffer and are characterized by the indicator 
specific dissociation constant K{ and the used indicator concentration i max . 



Adjustment to a specific experiment 

The idea of this semi-universal model for presynaptic calcium dynamics is to de- 
termine the universal parameters using single protein experimental data and to 
maintain the resulting values for the evaluation of different systems and experimen- 
tal conditions. Universality has to be understood as the postulated statement, that 
universal parameters are not the parameters that are most sensitive to a transfer 
from one experiment to another. The adjustment of the model to a specific exper- 
iment is achieved through the determination of the system-specific and condition- 
specific parameters. Note, that without any exception these are parameters with 
direct physiological interpretation. Therefore, most of them may be accessible up 
to a sufficient precision in several experiments. 

In the described procedure the main part of the model is determined by sources 
that are independent of the experiment under consideration. Therefore, the value 
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of the model is tested by its capability of reproducing the answer of the calcium 
concentration to single action potentials on the basis of those independent and fixed 
universal parameters. As in general not all specific parameters will have been de- 
termined in the experiment under consideration, we will fit the remaining unknown 
parameters to the single action potential calcium response. With the help of a 
thus defined model it should be possible to analyse the measured calcium transients 
evoked by LTD or LTP stimulation protocols. 

This has been executed for an experiment on single nerve terminals of pyramidal 
neurons in the neocortex of rats ||: Most specific parameters used in the model 
are directly accessible in this experiment. This concerns for example the form of 
the action potential applied to the nerve terminal, the geometry (i.e. the surface to 
volume ratio G = 3/pm), the concentration of endogenous buffer 6 max = 120/xM, 
and the characteristics and concentrations of the used calcium indicator (magnesium 
green). Therefore, we are in the situation that the whole model is a priori determined 
either by independent sources (concerning universal parameters) or by the available 
data from the experiment under consideration. The only unknown parameters are 
the surface densities of the calcium transport proteins. These are fitted to the 
measured calcium transient evoked by single action potentials (pu — 2.65/pm 2 , 
Pp = 4000//im 2 , and pe = Pp/30). The result Fig. |l| shows, that the measured 
calcium transient is reproduced correctly. 



80 



>- 20 

-20 
-40 
-60 
-80 

-50 50 100 150 200 250 

0.18 
0.16 
0.14 
0.12 
0.1 
0.08 

fa 

<] 0.06 
0.04 
0.02 
0.0 
-0.02 

-50 50 100 150 200 250 

t [ms] 

Figure 1: The calcium transients evoked by single action potentials in single boutons 
of pyramidal neurons in the rat neocortex (lOOpM magnesium green). The dotted 
line shows the experimental values (we thank H.J. Koester and B. Sakmann for 
kind permission to reproduce the data), and the full line shows the model result 
with fitted protein densities. 
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LTE stimulation 



In order to check if the thus defined model has predictive power, we calculate the 
intracellular calcium transients evoked by series of action potentials with varying 
frequency. Basically, the model parameters remain unchanged. As in the corre- 
sponding experimental setup (see || Fig. 9) a different action potential (compared 
to Fig. |l|) has been used, the action potential is adapted in the model and the channel 
densities are fitted to the single action potential calcium response, correspondingly 
(pu = 2.77/pm 2 , pp = 3400/pm 2 , and pe = pp/30). In addition the indicator 
concentration is increased, which is considerably higher in this experimental setup. 
Now the 10 Hz stimulus is applied and the result is shown in Fig. |2|. The model 




Figure 2: The calcium transients evoked by a 10Hz stimulus in single boutons of 
pyramidal neurons in the rat neocortex (500/zM magnesium green). The full line 
represents the model result and the dotted line the corresponding measured tran- 
sients || (we thank H.J. Koester and B. Sakmann for kind permission to reproduce 
the data). The dashed line represents the best single exponential fit to the experi- 
mental base line values IHl. 



result is in quantitative agreement with the calcium transients seen in the experi- 
ment: The intracellular calcium concentration reaches a new baseline level during 
the stimulation process which breaks down when the stimulus is switched off. The 
calcium concentration oscillates on the top of the new base line in coherence with 
the stimulation potential. The calcium reducing processes are slightly to slow in the 
model. 

On this basis we can calculate the calcium transients in response to stimulation 
protocols with various frequencies. We find that calcium transients do not overlap 
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for low frequencies typical for LTD-induction. The emergence of a new baseline 
in the calcium concentration at frequencies around 10 Hz may be interpreted as 
threshold for the induction of LTP. Note, that this threshold frequency strongly 
depends on the used calcium indicator concentration. This is especially relevant for 
the interpretation of experimental results. A stimulation with frequencies around 
50 Hz (typical for the induction of LTP) leads to a more pronounced enhancement 
of the calcium baseline. This qualitative behavior is in agreement with experiments 



carried out on dendritic spines of pyramidal neurons [26 



Conclusion 

Our new model for transients of the presynaptic intracellular calcium concentration 
evoked by various stimulation protocols reproduces the general behavior observed 
in experiment. It is exclusively constructed with parameters that have a direct 
physiological interpretation. Its basis are the single proteins properties. The char- 
acteristics of single proteins are considered to be universal in the sense that they 
remain unchanged for different experimental setups. The model has been adjusted to 
a specific experiment that measured intracellular calcium transients in nerve termi- 
nals of pyramidal neurons of the rat neocortex. To this end the parameters specific 
for this experiment have been extracted from it or, if not available, have been fitted 
to the single action potential calcium response. The model results turned out to be 
in quantitative agreement with the experiment. This applies not only to the presy- 
naptic calcium concentration response to single action potentials but also to 10 Hz 
stimuli. We did not find any reason for the involvement of calcium-induced-calcium- 
release in the induction of LTE. However, it seems that an additional mechanism 
(e.g. calcium channel inactivation) may be necessary to understand the induction of 
LTD on the level of presynaptic calcium transients. 

More generally, the separation of universal and specific parameters enables us 
to analyse different results observed in several experiments. With the help of the 
new model one may decide, if those differences are significant or due to different 
experimental setups. 
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